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ABSTRACT 



Context. Young supernova remnants (SNRs) exhibit narrow filaments of non-thermal X-ray emission whose widths can be limited 
either by electron energy losses or damping of the magnetic field. 

Aims. We want to investigate whether or not different models of these filaments can be observationally tested. 

Methods. Using observational parameters of four historical remnants, we calculate the filament profiles and compare the spectra of the 
filaments with those of the total non-thermal emission. For that purpose, we solve an one-dimensional stationary transport equation 
for the isotropic differential number density of the electrons. 

Results. We find that the difference between the spectra of filament and total non-thermal emission above 1 keV is more pronounced 
in the damping model than in the energy-loss model. 

Conclusions. A considerable damping of the magnetic field can result in an observable difference between the spectra of filament 
and total non-thermal emission, thus potentially permitting an observational discrimination between the energy-loss model and the 
damping model of the X-ray filaments. 
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1. Introduction 

Based on simple energetic considerations regarding the energy 
density of cosmic rays and the energy release per supernova ex- 
plosion, SNRs have long been thought to be sources of galactic 
cosmic rays. This presumption is supported by numerous detec- 
tions of non-thermal emission in the radio and X-ray band (e.g. , 
iKovama etalll995t ISlane et al.|[l999l l2001t iBamba et al.ll2000h 
observed from the direction of known SNRs and interpreted to 
be synchrotron radiation of relativistic electrons with energies 
uptolOOTeV. 

High-resolution observations of young SNRs performed 
with the Chandra satellite show that this emission of non- 
thermal radiation is concentrated in narrow regions on the limbs 
(IVink & Lamingll2003HBamba et al.l l2005). These regions of in- 
creased synchrotron emissivity close to the forward shock are 
called filaments and demonstrate the presence of high-energy 
electrons around their acceleration sites. 

The most plausible process for the acceleration of elec- 
trons is diffusive shock acceleration (P SA), which leads to a 
power-law distribution of particles (e.g.. Bell 1978; Blandford 
& Ostriker l 1978b iDrurvl Il983h . Although no clear evidence 
for relativistic-ion acceleration exists at shocks, the DSA- 
mechanism is also considered to be responsible for the accel- 
eration of cosmic -ray nuclei, as indicated by observations of 
non-relativistic ion acce leration at solar-wind shocks driven by 
coronal mass ejections (Rouillard et al.ll201ll) . However, many 
details of the DSA are still vague such as the maximum energy 
of particles, the role of the magnetic field, and how the particles 
are injected into the acceleration process (also referred to as the 
injection problem). Apparently, the investigation of the proper- 
ties of the non-thermal filaments may provide key information 
for a better understanding of the DSA-mechanism. In particu- 



lar, knowing the magnetic -field strength gives constraints on the 
maximum particle energy achievable in the acceleration process, 
which can help answering the question whether SNRs can ac- 
celerate particles to energies above the knee in the cosmic-ray 
spectrum. 

Accurate analysis of several SNRs shows that the filamentary 
structures are very thin compared with the radii of the remnants. 
This limitation of the filament widths is associated with a rapid 
decrease of the synchrotron emissivity that can be explained by 
energy losses of the electrons in a locally enhanced magnetic 
field. A number of authors have used that model to constrain the 
magnetic- field strength, the degree of turbulence and the obliq- 
uity (e.g.. iBamba et al]|2003t IParizot et all l2006t lArava et all 
2010). rArava et al.l d2010), for instance, have investigated several 
filaments of the remnant of Cas A and found that the magnetic 
fields of the filaments are highly turbulent and nearly perpendic- 
ular to the shock normal. Another important result of this and 
other studies is an estimate of the downstream magnetic-field 
strength that is higher than simple shock compression would 
suggest. Such observations indicate an additional amplification 
of the magnetic field in the shock region of the SNRs. A possi- 
ble amplification process could be a streaming instability i n the 
upstre am region as proposed by iLucek & Belli (2000) and iBelll 
(l2004t) . or the effects of preexist ing turbulent density fluctua tions 
on the propagating shock front (Giacal one & Jokipilll2007l) . 

Besides energy losses, also the magnetic field itself can 
limit the filament widths. Based on the turbulence relaxation 
downstream of the forward shock and neglecting any amplifi- 
cation process, iPohl et all (H005) have calculated that the tur- 
bulent magnetic field downstream can decay exponentially on a 
damping-length scale = 10 16 - 10 17 cm that is small enough to 
produce the narrow observable filaments. Furthermore, from ob- 
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servations of the post-shock steepening of the synchrotron spec- 
trum in Tycho's SNR it can be seen that also the damping of 
the m agnetic field fairly well desc ribes the corresponding X-ray 
data (Cassam- Chenai et al.l l2007). and thus, may appear within 
the filaments. 

Since the magnetic field controls the radiative cooling of 
electrons, high magnetic fields lead to strong cooling, and too 
few high-energy electrons remain capable of producing the 
gamma-ray emission, that is o bserved from regions n ear the 
edges of numerous SNRs (e .g.. lAharonian et al.ll2007k Acero 
et al. 120101 ; lAbdo et alJ^OllF Any gamma rays observed in 
such a case are likely to be hadronic in origin. Weak cool- 
ing leads to a large number of high-energy electrons and the 
possibility of gamma-ray emission through inverse Compton 
or bremsstrahlung processes. All these implications of the 
magnetic-field structure on the particle acceleration, gamma-ray 
emission and magnetic-field amplification make it necessary to 
understand the non-thermal filaments in detail. 

In this paper we investigate the properties of the filaments 
for both cases, filaments limited by electron energy losses or 
by damping of the magnetic field. For that purpose, using ob- 
servational values of some characteristic al SNR parameter, we 
calculate the X-ray emission of the filaments. The resulting fila- 
ment profiles then allow us to make specific predictions regard- 
ing the magnetic-field strength. We additionally calculate the to- 
tal non-thermal emission, which shall be referred to as "plateau", 
and compare their spectra with those of the filaments. It should 
be noted that in our models we only consider non-thermal syn- 
chrotron emission and restrict ourselves to the evolution of rel- 
ativistic electrons in the downstream region. Furthermore, we 
assume the electrons to be already accelerated at the shock front 
and treat our problem to be independent of the acceleration pro- 
cess. We also do not consider any electron propagation into the 
upstream region and simplify the SNRs to be spherical objects 
of constant downstream-velocity profile. Recent hydrodynami- 
cal simulations suggest that this oversimplification of a constant 
velocity is an acceptable approxim ation only for SNRs of a n age 
of less than several hundred years dTelezhinskv et al.ll2012h . im- 
plying that our models are restricted to SNRs being in the adia- 
batic phase and just entering the Sedov phase, respectively. 



where N(r, E) and P v (r, E) are the isotropic differential electron 
number density and the spectral emissivity per electron, respec- 
tively. Thus, we need to derive the appropriate electron distribu- 
tion within the filaments. 



2.1. The electron distribution 

In the following, we derive the electron distribution that is neces- 
sary to calculate the synchrotron emissivity. To do this, a trans- 
port equation describing the dynamics of a distribution of rel- 
ativistic electrons affected by advection, diffusion and energy 
lo sses needs to be solved. According to simulations described 
in iTelezhinskv et al.l d2012h . we can approximate the advection 
velocity of young SNRs to be constant downstream of the shock, 
implying that energy losses due to adiabatic deceleration can be 
neglected. Note that the non-thermal emission come from a thin 
spherical shell near the edges of the SNRs. If we restrict our 
treatment to a region near the shock that is crossed by the ad- 
vection flow on a timescale short compared with the age of the 
SNR, then we approximate the electron distribution with a one- 
dimensional steady-state solution. 

It is convenient to introduce a comoving spatial coordinate, 
Z — r s - r, where z = marks the position of the shock front 
at all times. Hence, restricting ourselves on the downstream re- 
gion, the one-dimensional transport equation for the isotropic 
differential number density, N = N{z, E), reads as follows: 



3_ 

dz 



oz 



^-^W(z,E)N] + Q(z,E) 
oz oE 



0. 



(3) 



In this equation v denotes the constant advection velocity of the 
electrons relative to the forward shock, D(z, E) is the diffusion 
coefficient, f3(z, E) - dE/dt represents the electron energy loss 
due to the emission of radiation, and Q(z, E) is the source term 
describing the injection of accelerated electrons into the propa- 
gation process. 

Since the electrons are likely accelerated by the DSA- 
mechanism at the forward shock (z - 0), we assume the elec- 
trons to be injected with a power-law dependence E~ s , where s 
is the injection index. Hence, the source term reads 



2. Modelling the filaments 

We calculate the X-ray intensity as a function of the projected 
radius, r p . It is an integral over the synchrotron emission coeffi- 

m j v dy. Using r = y+rt 
and taking into account that only emission originating inside the 
SNR contributes, the X-ray intensity can be written as 
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where r and r s denote the positions inside the SNR and the ra- 
dius of the SNR, respectively. Then, obtaining the corresponding 
spectrum of the filaments just involves an integration over the X- 
ray emission along the projected radius, whereas the spectrum of 
the plateau emission can be calculated as a volume integral over 
the synchrotron emission coefficient. In both cases, the isotropic 
synchrotron emission coefficient is needed, given by 



Mr) 



1 r 



N(r,E)P v (r,E)dE , 



(2) 



Q(z,E) = qoE s exp -— \6(z) , 

£cut/ 



(4) 



and includes a cut-off at energy E cu t, because the maximum pos- 
sible energy can be limit ed by either the fi nite accelera tion time 
of the SNR dDrurvll 19911) or energy losses dWebb et al.lll984l) . 

Eq. (0 can be solved using Green's method, implying that 
the solution can be written in terms of Green's function, 

dz' J dE' g(z,z';E,E')Q(z',E') 



') , (5) 
where Green's function g - g(z, z'\ E, E') satisfies 

-5{z-z')5{E-E') .(6) 
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Assuming the diffusion coefficient and the energy losses to 
be separable in a spatial and in an energetic part, 



D(z,E) = d(z)D(E) , 
p(z,E) = -b(z)B(E), 



(7) 
(8) 
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and that the spatial dependent terms are inversely proportional 
to each other, 



d(z)b(z) - a - const. , 

as well as introducing the substitutions 

G(z,z';E,E') 



B(E) 



(9) 



(10) 



g(z,z';E,E') = 

i ^ r dy 

Jo d(y) 

i r°° du 

a J E B(u) 

an analytical solution to Eq. (O can be found in the literature 
dLerche & Schlickeiserl [l980. see their Eq. (A20)) and reads 



(11) 



(12) 



G(x, x; A, A') 



@(A - A') 



2 <* "y JT, D(t) df 



x exp < - - 



[v(A - A') + x' - x] 2 
4j£/)(t)d* 



(13) 



where &(A - A') is the step function and A' = A(E'). Note, this 
analytical solution is valid only if Eq. ((9} applies. 

In general, the diffusion coefficient is an unknown param- 
eter. But it is often assumed that the diffusion proceeds in the 
so-called Bohm-regime. Hence, the diffusion coefficient can be 
written in the extreme-relativistic limit E » mc 1 as D — rir^c/3, 
where = E/(eB) and 77 > 1 are the gyroradius and gy- 
rofactor, respectively. Here, m is the mass of the electron, e 
is its charge, c is the speed of light, and B is the magnetic- 
field strength. Therefore, we may write D(E) = D^E, where 
Do = Tjc/(3eB). Additionally, we assume the emission of syn- 
chrotron radiation to be the main energy-loss process that is pro- 
portional to the square of the electron energy, B(E) = boE 1 , 
where bo = 4e 4 B 2 /(9m 4 c 7 ). Using these assumptions, Eq. (fT2l 
rewrites as 

1 f d« 1 

A(E)=- 

a J E 



so that 



and hence, 
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According to Eq. (1 1 Ob and Eq. ( 113b , Green's function then re- 
duces to 



g(x, x'\ E, E') = 



6(E' -E) 1 



2 ^InaboDoE 2 \In(f ) 



x exp 



SKA 



(14) 



Now, the electron distribution can be easily determined using 
Eq. (0 and the appropriate source distribution (0J, so that 



N(x(z) 



,E) = f dx' [ 
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Substituting n — E' /E, we can rewrite the desired result as 



N(x(z),E) 



/^OO — C 

r -^= 

Ji Vln(n) 



2 yfnaboDo J\ Vln(n) 

{ n£ f^(l-^)-4z)] 2 ) 
Xexp^-- : ^d«. (16) 



^cut 



2.2. Parameter used in the models 

The models we use are based on different assumptions on the 
magnetic field. In the first model, which shall be referred to as 
"energy-loss model", we assume the magnetic field to be spa- 
tially constant, whereas the second one, which is referred to as 
"magnetic-field damping model", includes the damping of mag- 
netic turbulence and assumes a spatially-dependent magnetic- 
field strength described by a profile following the relation 



B(z) = B min + (B 

n 



B 



rain) eX PI~^ 



(17) 



where is the damping length and z > 0. Here, we choose 
the minimum value of the magnetic field to be similar to that 
of the interstellar medium, B m [ n = 10 //G, whereas the maxi- 
mum value, B max , corresponds to the field at the shock. It should 
be noted that Eq. ([T71 > describes the averaged amplitude of the 
magnetic field in a given volume. Hence, we do not make any 
assumption on the magnetic-field direction and do not distin- 
guish between parallel and transverse diffusion, as is done in 
detailed calculations of diffusion coefficients (e.2., Marcowith 
et al. l2006h . 

According to the Rankine-Hugoniot conditions, the down- 
stream advection velocity can be expressed by the shock veloc- 
ity, v s , through v = v s /4, if we consider strong shocks with a high 
Mach number and a monatomic gas with adiabatic index of 5/3, 
leading to a shock compression ratio of 4. Thus, we neglect any 
non-linear effects expected to occur with efficient particle ac - 
celeration, which would modify the shock (Elli son et al.l l2004). 
Furthermore, the filament width, w, is defined as the length, at 
which the intensity described by Eq. ([TJ is reduced by a factor 
1 /e of its maximum. 

To have a realistic model, the values chosen for the shock 
velocity, filament width, and radius are based on reference val- 
ues of real SNRs. In our case, we consider the young remnants 
of the historical supernovae SN 1006, Cas A, Tycho and Kepler. 
It should be noted that the real filament widths found in the lit- 
erature have been measured in a certain X-ray band and not for 
an individual X-ray energy. However, as is shown in the calcu- 
lation done in Sect. [3] the shape of the filament profiles depends 
on the X-ray energy. Nevertheless, we relate the observational 
value of the width to a X-ray energy of 5 keV, since it is an en- 
ergy, at which no significant contribution from thermal emission 
is expected. 
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Additionally, we treat the cut-off energy to be the maximum 
electron energy that can be achieved in the acceleration process. 
Unlike cosmic-ray nuclei whose energy is probably age-limited 
due to the finite acceleration time available, we assume the maxi- 
mum electron energy to be loss-limited, since the electrons expe- 
rience synchrotron losses during their acceleration. By equating 
the acceleration time scale to the synchrotron loss time, it is pos- 
sible to derive an expression for the maximum electron energy 
in terms of the downstream magnetic-field strength and shock 
velocity (IParizot et al.ll2006h : 



£cut = £max - (8.3 TeV) J] 



-1/2 



( B(z = 0) 
\ 100 fiG 



-1/2 



1000 km/s 



■(18) 



Since diffe rent mechanism can account for magnetic-field am- 
plification ( iLucek & Bellll2000l : lGiacalone & Jokipiill2007h . the 
structure of the magnetic field is generally unknown within the 
shock region. We therefore simplify the problem by making the 
assumption that only the magnetic-field strength at the shock de- 
termines the maximum electron energy given by Eq. d 1 8b . In ad- 
dition, we have taken a shock compression ratio of 4. 

At last, we assume the injection index to be s = 2, which 
results from the Rankine-Hugoniot conditions for strong shocks 
dBellll978h . as well as Bohm diffusion (rj — 1), which implies the 
smallest possible value of the diffusion coefficient, as the mean 
free path of the particle is equal to the gyroradius. 

All parameter used are summarized in Table Q] Note, that it 
is possible that all four SNRs could exhibit similar shock veloc- 
ities. To take the uncertainties of this quantity into account, we 
perform the calculation of SN 1006 and Kepler for two different 
shock velocities. 



3. Results 

3. 1 . Energy-loss model 

In this model we treat the magnetic field to be constant, B(z) = 
B = const., implying no spatial dependence of the energy-loss 
term, b(z) = 1. Using d(z) = a/b(z), as well as Eq. (fTTT) . the 
spatial coordinate x(z) then scales as 

x(z(r)) = - = . 



Using the parameters given in Table Q~| we calculate the X- 
ray intensity as a function of the projected radius according to 
Eq. ([TJ. Here, the magnetic-field strength is a free parameter 
whose value can be chosen so that the filament widths match 
those found in the observations. In addition, we also determine 
the cut-off energy of the electron spectrum, which is connected 
to the magnetic field through Eq. d 1 8b . 

Reproducing the filament width for each SNR at a pho- 
ton energy of 5 keV determines the magnetic-field strengths as 
given in Table [2] For our examples the downstream magnetic- 
field strength ranges from about 100 yuG up to about 500 fiG. 
Remnants with narrower filaments exhibit a higher downstream 
magnetic field. These magnetic fields then imply cut-off ener- 
gies of the electron spectra in the energy range between 19 TeV 
and 37 TeV. The order of magnitude of the cut-off energies is in 
agreement with the results obtained from spectral modeling of 
the radio-to-X-ray spectra of young SNRs whose cut-off ener- 
gies of their electron distribution must be in the TeV-band, since 
the cut-off frequencies are generally found in the X-ray band 
dRevnolds & Keohandl 19991) . 

In Fig. (Q]) we illustrate the profiles of the filaments for 
four different photon energies calculated with the parameters of 



Tycho. To be noted from the figure is a frequency dependence of 
the filament widths, which can be explained by the energy loss of 
the electrons in a constant magnetic field and by the advection 
process. The advection length represents the distance covered 
by the electrons within the synchrotron loss time, T syn = E/\^\, 
and is given by 



'ad — VT syn — 



(2 x 10~ z pc) 



v s 9m c 
7 4e 4 B 2 E 
B 



200 fiG) \ 5000 km s" 



-)(2()!evf (19) 



Synchrotron radiation usually provides a continuum around a 
characteristic synchrotron frequency, 



3v L y 2 



3eBE 2 
47rm 3 c 5 



(1.3 x 10 18 Hz) 



B 



200 //G/ 1 20 TeV/ 



(20) 



where vl and y are the Larmor frequency and the Lorentz fac- 
tor of the accelerated electrons, respectively. Hence, the relation 
£ k y 1 ' 2 . Therefore, a dependence of the width on the frequency 
of the radiation of the form Z a( j oc y~'/ 2 would result. But the 
electrons are also affected by diffusion. With the Bohm diffusion 
coefficient, D = one thus obtains for the corresponding 

diffusion length 



dirt 



-4 



— ,(1.3xl0- P c) 



B 



,-3/2 



200 fxG ) 



,(21) 



which does not depend on the electron energy. Equating the ad- 
vection and diffusion length, the relation 



(31 TeV) 



(5OOO km s" 1 )! 



B 



200 fiG 



-1/2 



(22) 



can be derived. Electrons with energies E > E c can stream 
farther from the shock than advection alone would allow. 
According to Eq. (1201 and Eq.(l22b. the characteristic syn- 
chrotron frequency for the electrons with E > E c is then higher 
than 



VcCEc) - (3.1 x 10 1 



Hz)( 



5000 km s- 



(23) 



corresponding to a X-ray energy of about 13 keV, if v s = 
5000 km/s. Hence, the filament profiles at higher photon ener- 
gies, which need the most energetic electrons, show approxi- 
mately the same behaviour as can be seen in Fig. ((TJ for the 
5-keV and 10-keV profile. 

In addition, the advection and diffusion length also explain 
the relation between the filament widths and the corresponding 
magnetic fields given in Table [2] Narrower filaments require a 
faster decrease in the synchrotron emissivity, implying a shorter 
advection and diffusion length. And according to Eq. ( 1191 > and 
Eq. OTT) . this is given for higher magnetic-field strengths. 

Now, equipped with the X-ray intensity distribution estab- 
lishing the filament profiles, and the volume emissivity, we cal- 
culate the spectra of the filament and plateau for each of our 
examples. To obtain the filament spectrum, we integrate the in- 
tensity along the projected radius from r p = r s up to r p = r s - w. 
Because we do not know the electron source strength, qo, we 
are not interested in absolute fluxes. However, we can show 
the qualitative behaviour described by the appropriate photon 
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Table 1. Parameters used to derive the filament profiles and spectra. Furthermore, s = 2 and Bohm diffusion are assumed in all 
SNRs. 



SNR 


Age 

[yr] 


Distance 


r s 






v s 




w 




^min 


[kpc] 


Reference 


[arcmin]" 


[pc]* 


[km s" 1 ] 


Reference 


[arcsec] 


[pc]* 


Reference 


SN 1006" 


1000 


2.2 


1 


15 


10 


4900 


2 


20 


0.2 


3 


10 














(2900) 












Cas A 


330 e 


3.4 


4 


2.5 


2.5 


5200 


5 


1.5 


0.03 


6,7 


10 


Tycho 


440 


2.5 


8 


4 


3 


5000 


9 


5 


0.06 


6 


10 


Kepler-' 


410 


4.8 


10 


1.5 


2 


5040 


11 


3.5 


0.08 


6 


10 






(6.4) 






(3) 


(6720) 






(0.11) 







Notes. <D) Taken from lGreenl d2009). <6) Directly inferred using the distance and the appropriate quantity given in angular units. (c> Only used in 
the magnetic-field damping model. <rf) The shock velocity in the northwestern limb (value in brackets) is used, as well as in the no rtheastern limb, 
where strong electron accele ration appears to occur . (c) The supernova explosion may have been observed in 1680 ( Ashworth 1980). Otherwise, the 
detection of radioactive 44 Ti (Hartmann et al. 1997) and the analysis of the dynamics of this remnant (Fesen et al. 2006J) also suggest an explosion 
date in the late 17th century. Distance to this remnant is uncertain. Calculation is performed for the lower and upper limit (values in brackets) 
to the distance. 

Refe rences. (1) IWinkler et alj <l2003l) ; (2) iGhavamian et alj J2003) + iKatsuda et al.l (120091); (3) lBamba"et~ai] J2003I) : (4) lReed et afl < fl995T) ; 
(5) IVink et aO dl998h: (6)lBamba et all d2005h: (7)lArava et al.1 d2010l): (8)lTian & Leahvl d201 111: (9HHughesl d200Cl) + distance: ( 10) Revnoso & 
Goss d!999h : qi) IVinkU2008h + distance 



spectral indices, which should be sufficient for the comparison. 
Assuming the photon spectra to show a power-law characteris- 
tic, Af v = F v /hv cc y~ r , where F v and F are the flux and photon 
spectral index, respectively, we can then describe the spectra be- 
tween the photon energies hv and hV through 



r = 



(24) 



Here, we also want to calculate the differences between the pho- 
ton spectral index of the filament spectrum, Ff, and that of the 
plateau emission, F p , which are given at three different photon 
energies, E v , in Table [2] Additionally, in Fig. (O we show the 
photon spectral indices calculated with the parameters of Tycho. 

As can be seen from Table |2 as well as from Fig. (0, the 
spectra at higher photon energies rarely differ significantly. The 
plateau shows a steeper spectrum at lower photon energies. Up 
to a photon energy of 1 ke V the difference between the indices of 
filament and plateau is in the range 0.05 - 0.34, whereas at ener- 
gies higher than 1 keV the difference is always smaller than 0. 1 . 
This property can be explained by the effective radiation of the 
energetic electrons in the enhanced magnetic field. Accordingly, 
the most energetic electrons lose all of their energy inside the fil- 
aments, implying that the regions farther from the shock provide 
almost no contributions to the total emission of hard X-rays, so 
that both filament and plateau show nearly the same behaviour. 

In Table|2]we also show the indices of the filaments at three 
different photon energies. One can see that the parameters of 
Cas A, Tycho, Kepler and SN 1006 lead to the same spectral be- 
haviour, if their shock velocities are similar. However, using the 
shock velocities measured in the northwestern limb of SN 1006 
and resulting from the upper limit to the distance to Kepler, it 
turns out that the filament spectrum is softer and harder, respec- 
tively, than in the former case. 



3.2. Magnetic-field damping model 

As already mentioned above, in the model of magnetic-field 
damping we assume the magnetic-field strength to follow a pro- 



filament profiles in the energy-loss model 



CO 



0.1 keV 
1 keV 
5 keV 
10 keV 



1 ■ 




2.75 



Fig. 1. Non-thermal X-ray intensity as a function of the projected 
radius calculated for four different X-ray energies with the pa- 
rameters of Tycho given in Table [TJ as well as B — 310 fiG. The 
forward shock is located at r s = 3 pc. 



file described in Eq. (1 1 7b . which can also be written as 

Bmax ^min / Z 



Biz) = B n 



. -max "min 

1 + — — exp 



The spatial dependence of the energy-loss term then obeys the 
relation 



Hz) 



1 + 



fi eXP l 7 

D min \ 'd 



corresponding to a spatial variation of the diffusion coefficient 

-2 



diz) = a 



1 + 



■exp,-- 



because the product b(z)d(z) must be constant, as required by 
Eq. ||9). According to Eq. (11 Q , the spatial coordinate x(z) then 
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Table 2. Constraints on the downstream magnetic-field strength, cut-off energy of the electron spectrum, photon spectral index of 
the filament, and the difference between the photon indices of the spectra of filament and plateau at three different photon energies, 
E v , for four young SNRs. These values are calculated using the energy-loss model and the parameters given in TableQ] 



SNR 


B 
UiG] 


£cut 
[TeV] 




r f 






r P -r f 




E r = 0.1 keV 


E v = 1 keV 


E v = 10 keV 


E v = 0.1 keV 


E y = 1 keV 


Ey = 10 keV 


SN 1006" 


130 


36 


1.81 


2.33 


2.71 


0.31 


0.07 


0.02 




(110) 


(23) 


(1.99) 


(2.49) 


(2.94) 


(0.25) 


(0.05) 


(0.01) 


Cas A 


520 


19 


1.83 


2.33 


2.70 


0.31 


0.06 


0.01 


Tycho 


310 


24 


1.82 


2.33 


2.71 


0.32 


0.06 


0.01 


Kepler 6 


250 


26 


1.81 


2.32 


2.70 


0.31 


0.07 


0.01 




(230) 


(37) 


(1.75) 


(2.23) 


(2.61) 


(0.34) 


(0.08) 


(0.02) 



Notes. The values in brackets were calculated using the shock velocity v s = 2900 km/s as measured in the northwestern limb. ' ' The values in 
brackets were calculated using the upper limit of 6.4 kpc to the distance. 



spectral indices in the energy-loss model 




E v in keV 

Fig. 2. Photon spectral indices of the spectra of filament and 
plateau using the parameters of Tycho. 



scales as 



a f + 2 B 1 



^min) 



mm 



1 - ex P[-^ 



+2l d - 



B„ 



1 - exp 



■(25) 



Again, using the parameters given in Table Q~| we calculate 
the X-ray intensity as a function of the projected radius. In this 
case, the damping length and the maximum field strength, B max , 
are free parameters that can be chosen so that the filam ent widths 
match those observed. According to the calculation o f lPohl et al. 
(2005), the damping length should be in the range — 10 16 - 
10 17 cm (l& - 0.003 - 0.03 pc). To also investigate the influence 
of the damping length on the results, we perform the calculation 
using two different values of in each of our examples. Here, 
the larger value used for 1^ may describe weak magnetic-field 
damping, whereas the smaller one may cause a strong damping 
of the field. But note that energy losses are still included. 

Reproducing the filament width of each SNR at a photon 
energy of 5 keV requires the maximum field strengths given 
in Table [3] Depending on the damping length, the maximum 
field strength can be found in the range between 50 fiG and 
260 fiG, implying, according to Eq. ( fTSI ). cut-off energies be- 
tween 27 TeV and 62 TeV. Furthermore, it turns out that an in- 
creased damping length requires an increased field strength, be- 



cause the electrons radiate efficiently in a larger volume, result- 
ing in wider filaments. To retain the observed filament widths, 
it is then necessary to have a higher magnetic field that, on the 
other hand, also leads to a smaller cut-off energy of the electron 
distribution. 

If the damping length is too small, the observed filament 
widths cannot be realized for any maximum field strength. In 
these cases the intensity first decreases but then increases again 
even for the 5-keV profile, so that the typical shape of the fila- 
ment profiles is not given anymore. Therefore, we use damping 
lengths in the case of strong damping, for which the profiles just 
still exhibit the typical shape. For instance, using the given pa- 
rameters of SN 1006, the damping length used in the calculation 
should not be smaller than 0.02 pc. 

The filament profiles calculated with the parameters of 
Tycho for four different photon energies are illustrated in Fig. 
(0. As can be seen from the figure, there is also a dependence of 
the filaments on the frequency of the X-rays. This dependence 
is based on the spatial variation of the magnetic -field strength. 
Only in regions very close to the shock front the electrons radi- 
ate in fields of high magnitude, so that even the most energetic 
of them can emit photons of several keV in energy only in a 
small volume. The electrons remain energetic when they prop- 
agate into the downstream region, where they radiate at lower 
frequencies. Hence, we expect increased emission of low-energy 
X-rays in regions farther from the shock. For instance, using the 
cut-off energy of the electron spectrum in Tycho obtained for the 
case of strong damping, E cu i = 45 TeV, one finds that accord- 
ing to Eq. d20t . even the most energetic electrons located in a 
magnetic field of B — 10 fiG have their synchrotron continuum 
around the characteristic frequency v c (E cu {) = 3.3 x 10 17 Hz, 
corresponding to about 1.4 keV in X-ray energy. This issue can 
be seen for the 0.1 -keV and 1-keV profile in Fig. (01. The X-ray 
intensity does not decrease with decreasing projected radius as 
happens in the energy-loss model, but remains nearly constant 
and even increases, respectively. 

Using the X-ray intensity distribution, as well as the volume 
emissivity, we now calculate the spectra of filament and plateau 
for each example. Again, we integrate the intensity along the 
projected radius between r p = r s and r p = r s - w in order to 
obtain the filament spectrum. The difference between the photon 
spectral indices of filament and plateau, Ff - Fp, at three photon 
energies is given in Table [3] whereas in Fig. (0| we illustrate the 
photon spectral indices calculated with the parameters of Tycho. 

As can be seen again, the spectrum of the plateau is steeper 
than that of the filament. However, the difference between the 
spectra of filament and plateau depends on the chosen damp- 
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0.1 keV 
1 keV 




t , , , , U 

2.75 2.8 2.85 2.9 2.95 3 

r p in pc 



Fig. 3. Non-thermal X-ray intensity as a function of the projected 
radius calculated for four different X-ray energies with the pa- 
rameters of Tycho given in TableQ] as well as = 0.008 pc and 
Bmax = 85 yuG. The forward shock is located at r s = 3 pc. 

ing length. At relatively large damping lengths (weak damping) 
the difference between the indices over the hole energy range 
is smaller than 0.1, whereas at smaller damping lengths (strong 
damping) it also takes values in the range 0.1-0.2. The small dif- 
ferences at larger damping lengths are due to the higher magnetic 
fields that need to be chosen in order to retain the observed fila- 
ment widths. Hence, considerable energy losses have to be taken 
into account. Similarly to the energy-loss model, this results in 
an almost equal behaviour of the spectra of filament and plateau 
at high photon energies. In contrast, the magnetic fields used at 
smaller damping lengths are low enough to result in spectra that 
show significant differences among each other. 

Finally, one can also see from Table[3]that the filament spec- 
trum becomes steeper with decreasing damping length, in partic- 
ular at small X-ray energies. This is due to the lower magnetic- 
field strengths used in that case. Although the weaker magnetic 
fields imply higher cut-off energies, which harden the spectra, 
their influence is not sufficient enough to result in spectra simi- 
lar to those found at larger damping lengths. Besides, as in the 
energy-loss model and independently of the damping length, Cas 
A, Tycho, Kepler and SN 1006 show roughly the same spectral 
behaviour, if the shock velocities are similar, whereas the fila- 
ment spectrum obtained from the shock velocity of the north- 
western limb of SN 1006 has a steeper profile. In contrast, the 
upper limit to the distance to Kepler implies a harder filament 
spectrum. 

4. Conclusions 

Compared to the magnetic-field damping model, the spectra of 
filament and plateau obtained in the energy-loss model exhibit 
larger spectral indices. This can be explained by the considerable 
energy losses leading to the evolution towards a softer electron 
distribution in the energy-loss model, and hence, resulting in X- 
ray spectra that are softer than those obtained in the damping 
model. 

In case of a weak magnetic-field damping the difference be- 
tween the spectral indices of filament and plateau over the full 
X-ray spectrum is smaller than 0.1, which is probably to small to 
be detectable. Only if there is a strong damping, our calculation 



filaments in young supernova remnants 

spectral indices in the damping model 



X3 




E v in keV 

Fig. 4. Photon spectral indices of the spectra of filament and 
plateau using the parameters of Tycho. 



suggests a measurable difference between the spectra in some 
SNRs, since the difference between the spectral indices of fila- 
ment and plateau can take values of almost 0.2 at X-ray energies 
higher than 1 keV. 

In the energy-loss model the difference between the indices 
of filament and plateau above the X-ray energy of 1 keV is even 
smaller than 0.1, so that a possible detection can be excluded 
here, too. On the other hand, the difference between the indices 
below 1 keV is larger than 0.1, and at a photon energy of 0. 1 ke V 
it is even approximately 0.3. This might suggest that there is 
a measurable difference in the spectra of filament and plateau 
at small X-ray energies, if the filaments are limited by energy 
losses. However, on account of the interstellar photoelectric ab- 
sorption of the soft X-rays, these different spectral characteris- 
tics are probably not detectable, too. Furthermore, the plasma 
downstream of the forward shock is at high temperature, imply- 
ing also thermal emission contributing to the soft X-ray band, 
and thus, complicating a clear identification of the non-thermal 
emission. 

Hence, if there is no measurable difference between the spec- 
tra of filament and plateau, it is not possible to make definite 
predictions from the comparison of the spectra whether the fil- 
aments are limited by energy losses of the radiating electrons 
or by damping of the magnetic field. But if a significant differ- 
ence appears, our calculations then suggest that the filaments are 
limited by the magnetic field itself. 

It should be noted that our results presented here have been 
derived using Bohm diffusion. According to Eq. (1211 i. a larger 
diffusion coefficient with gyrofactor rj > 1 would imply a larger 
diffusion length, resulting in significant widening of the fila- 
ments, because now, the regions farther from the shock contain 
a sufficient number of high-energy electrons contributing to the 
intensity. Widening must then be compensated by a higher mag- 
netic field to retain the observed filament widths. Hence, the 
magnetic-field strengths derived in our models represent lower 
limits for the chosen parameters. The calculation then shows 
that a larger diffusion coefficient results in softer spectra due 
to a lower cut-off energy, which decreases with increasing gy- 
rofactor. However, the final results regarding the differences in 
spectral indices do not change fundamentally. 

In a last step we want to compare the predictions derived 
here with observations. At first, we notice that, independently of 
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Table 3. Constraints on the maximum magnetic-field strength, cut-off energy of the electron spectrum, photon spectral index of the 
filament, and the difference between the photon indices of the spectra of filament and plateau at three different photon energies, E v , 
for four young SNRs. These values are calculated using the magnetic-field damping model and the parameters given in Table [TJ 



Weak Damping 



SNR 


h 


Bmax 


E C ut 






If 








r f -r p 








[pc] 


IpG] 


[TeV] 


Ey 


= 0.1 keV 


E v = 1 keV Ey 


= 10 keV 


Ey 


= 0.1 keV 


E v = 1 keV 


Ey 


= 10 keV 


SN 1006 


0.03 


65 


50 




1.78 


2.12 


2.64 




0.04 


0.08 




0.07 




(0.03) 


(57) 


(32) 




(1-93) 


(2.35) 


(2.95) 




(0.05) 


(0.09) 




(0.08) 


Cas A 


0.015 


260 


27 




1.71 


1.98 


2.51 




0.04 


0.02 




0.01 


Tycho 


0.02 


150 


34 




1.71 


2.00 


2.57 




0.06 


0.04 




0.02 


Kepler 


0.03 


135 


36 




1.71 


2.01 


2.59 




0.05 


0.04 




0.02 




(0.03) 


(115) 


(52) 




(1.67) 


(1.93) 


(2.43) 




(0.04) 


(0.05) 




(0.02) 


Strong Damping 


SNR 


k 


Smax 


Ecut 






r f 








r f -r p 








[pc] 


IpG] 


[TeV] 


Ey 


= 0.1 keV 


Ey = 1 keV Ey 


= 10 keV 


Ey 


= 0.1 keV 


E v = 1 keV 


Ey 


= 10 keV 


SN 1006" 


0.02 


64 


51 




1.81 


2.16 


2.66 




0.04 


0.10 




0.13 




(0.02) 


(56) 


(32) 




(1.97) 


(2.39) 


(2.96) 




(0.04) 


(0.11) 




(0.13) 


Cas A 


0.004 


115 


40 




1.81 


2.12 


2.58 




0.07 


0.16 




0.19 


Tycho 


0.008 


85 


45 




1.82 


2.14 


2.61 




0.05 


0.13 




0.17 


Kepler 6 


0.01 


80 


47 




1.81 


2.14 


2.62 




0.04 


0.11 




0.14 




(0.012) 


(80) 


(62) 




(1.74) 


(2.03) 


(2.48) 




(0.03) 


(0.09) 




(0.12) 



Notes. (fl) The values in brackets were calculated using the shock velocity v s = 2900 km/s as measured in the northwestern limb. {b) The values in 
brackets were calculated using the upper limit of 6.4 kpc to the distance. 



the model, the parameters from the remnants of Cas A, Tycho, 
Kepler and SN 1006 lead to nearly the same spectral behaviour 
in case of similar shock velocities, as can be seen from the spec- 
tral indices in Table [2] and Table [3] However, th e analysis of 
the filamen t spectra of these remnants reported bv lBamba et alJ 
J200ll2005h reveals significant differences among the spectral 
indices obtained from the fit of an absorbed power-law model. 
Compared to our spectra whose calculation has been done using 
an injection index resulting from an unmodified shock (s = 2), 
the observation may be an indication for different electron in- 
jection indices in these remnants, implying shocks that are dif- 
ferently affected by non-linear effects due to differences in effi- 
ciency in the particle acceleration. 

Regarding the magnetic-field strengths, we take, as an exam- 
ple for comparison with our results, th e non-thermal filaments of 
Cas A analysed bv lArava et al.l (J2010). From the best-fit param- 
eters used to fit the observed filament spectra, the magnetic field 
has been derived to be in the range (30 - 70) fiG. These values 
are consistent with those derived from the magnetic-field damp- 
ing model, in which the magnetic field varies, according to Eq. 
( 117l i and the values from Table [3] between the field strengths 
(10 - 260) yuG for weak damping and (10-1 15) fiG for strong 
damping, respectively. For comparison, the constant magnetic 
field derived from the energy-loss model is several times higher, 
B = 520 fiG. This might suggest that the non-thermal filaments 
of Cas A are limited by the damping of the magnetic field. 

Another comparison concerns the magnetic fields in SN 
1006 and T ycho. Using the dat a from radio up to TeV- 
observations, lAcero et al.l (1201 Oh have analysed the multi- 
wavelength spectrum of SN 1006 in the framework of a lep- 
tonic and hadronic origin for the gamma-ray emission, giving a 
magnetic-field of ~ 30 fiG in the leptonic model and a magnetic 
field of ~ 120 fiG in the hadronic model, respectively. Moreover, 
combining radio and X-ray data with recent TeV-observations 
performed with the VERITAS instrument, the magnetic field of 
Tycho has been estimated to be ~ 80 fiG in a leptonic-dominated 
model, whereas a hadronic dominated model yields a magnetic 
field of ~ 230 fiG (lAcciari et al.ll201 H) . Compared to our model 



predictions given in Table [2] and Table [3] we notice that the 
magnetic fields derived from the energy-loss model are in good 
agreement with those estimated from the hadronic model used to 
describe the observed spectra of SN 1006 and Tycho. In contrast, 
the predictions from the magnetic-field damping model suggest 
the leptonic model for the origin of the gamma-ray emission 
from these remnants. It should be noted that current gamma-ray 
observations do not reach the spatial resolution of those done in 
X-rays, so that the magnetic fields estimated using gamma-ray 
observations of SN 1006 and Tycho are averages over a region 
much larger than the filaments, implying that the observed val- 
ues do not necessarily match those found for the filaments. 

To discriminate between the energy-loss model and 
magnetic-field damping model, and hence between a leptonic 
and a hadronic origin of TeV-band gamma-ray emission, one 
may either search for differences between X-ray spectra of fil- 
aments and plateau, as calculated in this paper, or perform 
gamma-ray observations with higher spatial resolution. 

Acknowledgement 

We acknowledge support by the "Helmholtz Alliance for 
Astroparticle Phyics HAP" funded by the Initiative and 
Networking Fund of the Helmholtz Association. 

References 

Abdo, A. A., Ackermann, M., Ajello, M., et al. 201 1, ApJ, 734, 28 
Acciari, V. A., Aliu, E., Arlen, T., et al. 2011, ApJ, 730, L20 
Acero, R, Aharonian, R, Akhperjanian, A. G., et al. 2010, A&A, 516, A62 
Aharonian, R, Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2007, ApJ, 661, 
236 

Araya, M., Lomiashvili, D., Chang, C, Lyutikov, M., & Cui, W. 2010, ApJ, 714, 
396 

Ashworth, Jr., W. B. 1980, Journal for the History of Astronomy, 11, 1 
Bamba, A., Koyama, K., & Tomida, H. 2000, PASJ, 52, 1157 
Bamba, A., Yamazaki, R., Ueno, M., & Koyama, K. 2003, ApJ, 589, 827 
Bamba, A., Yamazaki, R., Yoshida, T., Terasawa, T., & Koyama, K. 2005, ApJ, 
621, 793 

Bell, A. R. 1978, MNRAS, 182, 147 



8 



R. Rettig & M. Pohl: Non-thermal X-ray filaments in young supernova remnants 



Bell, A. R. 2004, MNRAS, 353, 550 

Blandford, R. D. & Ostriker, J. P. 1978, ApJ, 221, L29 

Cassam-Chenai, G., Hughes, J. P., Ballet, J., & Decourchelle, A. 2007, ApJ, 665, 
315 

Drury, L. O. 1983, Reports on Progress in Physics, 46, 973 
Drury, L. O. 1991, MNRAS, 251, 340 

Ellison, D. C, Decourchelle, A., & Ballet, J. 2004, A&A, 413, 189 
Fesen, R. A., Hammell, M. C, Morse, J., et al. 2006, ApJ, 645, 283 
Ghavamian, P., Winkler, P. R, Raymond, J. C, & Long, K. S. 2002, ApJ, 572, 
888 

Giacalone, J. & Jokipii, J. R. 2007, ApJ, 663, L41 

Green, D. A. 2009, Bulletin of the Astronomical Society of India, 37, 45 
Hartmann, D. H., Predehl, P., Greiner, J., et al. 1997, Nuclear Physics A, 621, 83 
Hughes, J. P. 2000, ApJ, 545, L53 

Katsuda, S., Petre, R., Long, K. S., et al. 2009, ApJ, 692, L105 
Koyama, K., Petre, R., Gotthelf, E. V., et al. 1995, Nature, 378, 255 
Lerche, I. & Schlickeiser, R. 1980, ApJ, 239, 1089 
Lucek, S. G. & Bell, A. R. 2000, MNRAS, 314, 65 
Marcowith, A., Lemoine, M., & Pelletier, G. 2006, A&A, 453, 193 
Parizot, E., Marcowith, A., Ballet, J., & Gallant, Y. A. 2006, A&A, 453, 387 
Pohl, M., Yan, H., & Lazarian, A. 2005, ApJ, 626, L101 
Reed, J. E., Hester, J. J., Fabian, A. C, & Winkler, P. F. 1995, ApJ, 440, 706 
Reynolds, S. P. & Keohane, J. W. 1999, ApJ, 525, 368 
Reynoso, E. M. & Goss, W. M. 1999, AJ, 118, 926 
Rouillard, A. P., Odstrcil, D., Sheeley, N. R., et al. 2011, ApJ, 735, 7 
Slane, P., Gaensler, B. M., Dame, T. M., et al. 1999, ApJ, 525, 357 
Slane, P., Hughes, J. P., Edgar, R. J., et al. 2001, ApJ, 548, 814 
Telezhinsky, I., Dwarkadas, V. V., & Pohl, M. 2012, Astroparticle Physics, 35, 
300 

Tian, W. W. & Leahy, D. A. 201 1, ApJ, 729, L15 
Vink, J. 2008, ApJ, 689, 231 

Vink, J., Bloemen, H, Kaastra, J. S., & Bleeker, J. A. M. 1998, A&A, 339, 201 

Vink, J. & Laming, J. M. 2003, ApJ, 584, 758 

Webb, G. M., Drury, L. O., & Biermann, P. 1984, A&A, 137, 185 

Winkler, P. E, Gupta, G, & Long, K. S. 2003, ApJ, 585, 324 



9 



